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Abstract 

We provide the mathematical structure and a simple, transparent and 
rigorous derivation of the magnons as elementary quasi-particle exci- 
tations at low temperatures and in the infinite spin limit for a large 
class of Heisenberg ferromagnets. The magnon canonical variables 
are obtained as fluctuation operators in the infinite spin limit. Their 
quantum character is governed by the size of the magnetization. 
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1 Introduction 



The appearance of spin waves in quantum ferromagnets at low temperatures 
is one of the most basic physical quantum characteristics of quantum spin 
systems. It amounts to a boson representation of the low temperature ele- 
mentary excitations of a spin system. The basic steps in the understanding 
of this phenomenon are done by Bloch QXJ , Holstein - Primakoff [[|, van Kra- 
mendonk - van Vleck H, and in the more technical work by Dyson [f|,|5|]. 

From the point of view of mathematical physics, one discovers rigorously spin 
wave properties at regular times as upper or lower bounds of correlations 
for low temperature or ground states (see e.g. |J). Undoubtly, the so-called 
Bethe ansatz |7J] is the most representative low temperature model of the spin 
wave theory, and, as is well known, to prove or disprove that the Bethe ansatz 
is correct for some models is a serious uptodate activity in mathematical 
physics. On the other hand, for a long time the spin wave theory of Holstein 
- Primakoff itself calls for a simple, transparent and mathematically rigorous 
setting, in the sense that the conditions for the appearance of spin waves are 
clearly formulated, and that then the derivation of the spin waves or magnons 
is rigorously obtained. 

As far as we know, the most serious attempt for this has been made in ||. 
On the basis of a classical domination principle which makes clear how to 
create the situation of a unique ferromagnetic quantum ground state, these 
authors define a 'physical' Hamiltonian, which is of direct relevance to the 
analysis of the Holstein - Primakoff and the Dyson formalism. They construct 
a 'quadratic boson Hamiltonian' in configuration space, and modulo an ap- 
proximation, called plausible for large spins, they are able to obtain satisfying 
results on the asymptotic equality of the free energies of the quadratic bo- 
son Hamiltonian and the Heisenberg model in any dimension for low enough 
temperatures. 

Here, stimulated by , we are able to clarify a number of aspects of the spin 
wave or magnon theory which were still absent. As in M, we limit ourself 
to ferromagnets. Our main results are that we are able to use mathematical 
results on non-commutative central limit theorems in order to scrutinize the 
large spin limit correctly and to give a rigorous scheme for the formation of 
bosons. We are able to perform this programme without any uncontrollable 
approximation. The result is that the magnon canonical variables are noth- 
ing but fluctuation operators, whose mathematical structure is developed 
in ||S , TTJ] . We perform rigorously the infinite spin limit, and prove that the 



Heisenberg model at low temperature becomes a system of non-interacting 
magnons. We discuss the quantum character of the magnons as a function 
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of the magnetization at low temperatures. 

Our results are on the level of the equilibrium states going beyond the study 
of its thermodynamic properties. We create the right conditions in order 
that the Heisenberg model system converges to a system of magnons. All our 
conditions are in the direction of enough ferromagnetism. We are convinced 
that results can be derived along the lines of what we have here also for 
antiferromagnets if one is able to pin the suitable conditions on the interaction 
constants in order to get specific antiferromagnetism. 



2 Low temperature magnons 

At each site x of a finite domain A of a cubic lattice Z", consider the spin- 
(2S + 1) variables S 1 (x), S 2 (x), S 3 (x), in the representation on <&f = _ s (C 2 )^ 
given by 

S"(x) = £ o?(s), (1) 

t=-S 

H = 1,2,3, the crf(x) Pauli matrices satisfying 

[al{x),a]{y)} = 2i5 id 5 x , v af(x), x,yel? (2) 
and its cyclic permutations of the components (1, 2, 3). Let 

±, , _ crl(x)±ia?(x) 

°i \ X ) — 2 

and 

S±(x) = £ ot(x). 

i=-S 

The Heisenberg model Hamiltonian on A is given by 

H A = -~ E {j(x,y)[S 1 (x)S 1 (y) + S 2 (x)S\y)] + J 3 (x,y)S 3 (x)S 3 (y)} 

x,yeA 

+hJ2S 3 (x). (3) 

We assume finite range, translation invariant interactions, i.e. J(x, y) = 
J(\x — y\), J 3 (x,y) = Js(\x — y\) and J(z) = J^{z) = for z large enough, 
we assume also J(x,x) = J 3 (x,x) = 0. 

Remark that the representation of the spin variables is special in the sense 
that the S^(x) ([!]), i.e. per lattice point x, are permutation invariant for 
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arbitrary permutations of the spin index i = —5, —5 + 1, . . . , 5. The Hamil- 
tonian on the other hand is not permutation invariant under permutations 
of the lattice indices igZ". 

Conform to Bloch, Holstein - Primakoff , Dyson , we are interested in the 
5 tending to infinity limit. In order to keep the model (^) thermo dynamically 
stable, one has to rescale it by the factor 25 + 1: Ha —> iff = (25 + 1) _1 //a- 
Applying a rescaled magnetic field h — > (25 + l)h, we get 

Hi = -W^—T\ E {j(x,y)[S 1 (x)S 1 (y)+S 2 (x)S 2 (y)] 

+J 3 (x,y)S 3 (x)S 3 (y)} + hY,S 3 (x). 

Or, rewritten, using J(x, y) = J(y, x) and J(x, x) = 0: 

= E {2J(x,y)[S + (x)S-(y) + S-(x)S + (y)] 

+J 3 (x,y)S 3 (x)S 3 (y)} + hY,S 3 (x) 

x£A 

= - 2(2 o + 1 ) E {4J(x,y)S + (x)S~(y) + J 3 (x,y)S 3 (x)S 3 (y)} 

+hY,S\x). (4) 

xeA 

We are interested in the equilibrium states of this model in the 5 — > oo 
limit, and in the thermodynamic limit A — > 7L V . For each finite volume A, 
iff is permutation invariant for one-site observables a\ (x) in spin space. Let 
uj be an equilibrium state, then per lattice site x 6 A, the state uj is again 
permutation invariant. In order to be able to give a definite mathematical 
meaning to limits of the type, for one-site x G A, 

lim^L, 

s-+oo V25 + 1 

in the sense of non-commutative central limits, the state uo has to satisfy a 
cluster property 

lim uj (AAx)BAx)) = uj(A)uj(B), 

i—>oo 

where A(x) and B{x) are products of the Pauli matrices (see e.g. ||). How- 
ever all clustering permutation invariant states are product states. This im- 
plies that uj, the infinite volume, infinite 5 limit equilibrium state at inverse 
temperature (3 of iff, is a product state in spin space: 

uj{AiBj) = uj{A)uj(B) 
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for i ytz j and A, B one-site observables in spin space, i.e. observables on 



B^aCC 2 )*. 



First we remark that for this product state in spin space, one has that in the 
o;-weak topology |§: 

weak - Jim i^ 2 g + 1 E a l( x ) \ = w(<r\x)) = ^(a tt ). 

The last equality follows from space translation invariance, which is a conse- 
quence of the space translation invariance of (^|). 

For k G A* = {k = ?fn ; n G Z^}, |A| = Z7, let: 



±ik.x 



ik.x 



and 
Then 



I I x-eA 

^(*) = 7T^E(^)- W (°?( a; ))) e ' 
I I x-eA 

dt(k) = a>(k) + \\\ l >M"')Sk,o- 



are fluctuation operators in spin space and in volume space. The infinite S 
limit of these operators is known to exist due to the product state character 
of oj @, i.e. the following limits exist: for all A G R, 

lim uj (exp{iX\FUk) + FUkY 

S—*oo V 

and 



lim u (exp{X[Fl(k) - Fl(k)*}}) . 



The limit operators lim5_ >00 Fg(k) are denoted F$(k), they still depend on 
the volume A. It is straightforward to check that these limits satisfy the 
canonical commutation relations. 

Before taking the infinite iS-limit, the operators 
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can be used to rewrite iff (^), after Fourier transformation, in the form 

+/ i (|A|(25+l)) 1 / 2 ^|(0), 

where J(Jfe) = £ z J(z, 0)e" ifc - z , J 3 (Jfe) = £ 2 J 3 (z, 0)e" ifc - 2 . 

The following Lemma characterizes the rotational invariance around the third 
axis, if there is enough ferromagnetism. 

Lemma 1 For h sufficiently large, i.e. for 

h > £[J 3 (x,0) - J (^°)] > °! 
x-eA 

any equilibrium state uj satisfies in the A — > Z u , S — > oo /zmii: 

a;((T ± ) = 0. 

Proof. Compute 

= (2STT)M EE{w(*,»)«fW«5-M 



i=-S 1,3/ 

aT 



-J 3 (x,y)fo (i)trj(y) + ajfo)^ (s)]} - —J^. ( x ) 



Then the time invariance of u implies 

= lim \wau([Hl,±rJ2^( x )i) 

= - lim m £{2J(z,i/)a;(0- 3 )u;(O -2J3{x,y)u{a 3 )u(a-)\ 

|A|-»oo |A| X J/ >• > 

-2/ia)(<r"), 



or 



a") jcu(a 3 ) £[J 3 (*, 0) - J(*, 0)] -h\ = 0. 
Since —1 < uj(a 3 ) < 1, taking ft, > 2,JJ 3 (z, 0) ~ J( z ify] > ensures 

^l^oi-^ou-ko. 
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Hence 

cj(cr") = 0. 



□ 

This result means also that the operators Fg and J-~g, as defined in (|5|) and 
@ coincide: 

Tf{k) = Ff(k), 
and hence also in the infinite S'-limit: 

jF±(jfe) = F ± (fc). 
We compute the commutators for finite S: 



|A| 1 /2( 2l S + l)i/2^ (/c g) ' 



and: 



2 

= ± |A|V2(25+1)V2^ (9±A:) - 
Therefore the limits satisfy the boson commutation relations: 

\F + {k)^-{q)] = lim = ^(a 3 )5 M , (7) 

and 

[^ 3 (A;),^(g)] = lim [F%{k),jf(q)] = ±2u{a ± )5 Kq = 0, (8) 

S— >oo 

on the basis of Lemma |l|. 

The J r± {k) are the above rigorously defined magnon creation and annihila- 
tion operators. Their existence and explicit properties are established as a 
straightforward application of the non-commutative central limit theorems 
in [|9] , and the condition of ferromagnetism in Lemma [I]. 

Now we proceed by determining the equilibrium state u completely. We use 
the definition of equilibrium state by means of correlation inequalities ITT 
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i.e. for temperatures T > {(3 = < oo), a state u> is an equilibrium state 
if and only if it satisfies the energy-entropy balance inequalities: 

lim lim fiu, (X*[Hl,X\) > oo(X*X)ln^^-, (9) 

for all local observables X. We prove: 

Theorem 1 In the ferromagnetic region (see Lemma\^), in the infinite S- 
limit, the equilibrium state uj of the Heisenberg model is a quasi-free state 
on the fluctuation operators algebra, generated by the {J^(k),k G [0,271-]^}, 
completely determined by the two-point function 



(f + (q)F-(q)) = pm ^ ){J ^J m+h] _ ■ (10) 



Proof. Take first X = JF 5 (g), and compute 
= -Y,2J(km(k),F s (q)]F s (k) 

k 

-Y,\Uk) ([n(k),^s(Q)]^l(-k)+n(k)[n(-k)^s(Q)] 
k z 

+h\A\^(2S + l) 1 / 2 [Fl(0),Fs(q)} 

+ |A|i/2( 2 g + l)V2 E Uk) (Fsi<l ~ k)Fl(-k) + Fl(k)F s (Q + k 
-2hTs{q). 



Then: 



~ WWTW 2 ? J{k)u ^ {q):Fl{k ~ q)Ts 



|A|i/2( 2 g+ 1)V2 E MkP {Tt{q)Ts{q - k)Fl(-k) 

|A|i/2( 2 g + 1)V2 E MQu {Tt{q)Tl{k)Ts{q + k) 
2(3hu (Tt{q)Ts{q)). 



Use that for A, B one-site observables in spin space such that 00(A) = 
(see !): 

(2S + ^ ^{AYTsiB^siA)) = u{A*A)u{B), 

to calculate 

Ymvju (Tt{q)[HtTs{q)]) 
lA L E J(k) Jim a; (j${q)Fj(k)) u(a 3 )6 k , q \A\^ 



5^oo 

2(3 



+ ^ E J 3(A:) jim w (J$(q)rj(q + kj) a;(a 3 )4,o|A| 1/2 



A- 



2(3h lim w fe(g)^ s 

6— »oo V 



Jim a; (^ s + (g)^ s (g)) = a; (^ + (g)^-(g)) . 



and 
Then 

lim lim (3uo(Ft{q)[HtFs{q)\ 

A— >ooo— >oo v 

= 2(3 {o;(a 3 )[J 3 (0) - J(q)} - h) u (^ + (g)^(g) 
After substitution into the correlation inequality (||), one gets: 

Interchanging the role of J r+ {q) and JF~(g), i.e. take now X = J 7 s(q) in ([ 
and repeat the computation above, then: 

2 , { -,,, )Wc) - JWlrt} > la WMj 

These two inequalities, combined with the commutation relation (0) yield 
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or the expected two-point function 



e 2/3{-^( CT 3)[J 3 (0)-J( g )]+fe} _ I ■ 



Finally if one takes for X higher order monomials in the J~s(q), one derives 
readily also from (0), that the higher order point correlation functions are 
sums of products of this two-point function, proving that the state u is quasi- 
free. As this amounts to a straightforward computation, we leave it as an 
exercise for the reader. 

□ 

The basic two-point function ( |i0|) contains still the magnetization u(o~ 3 ). 
Using (H) and @, one gets a self-consistency equation for the magnetization 

* £ w (^(,)^ ta )) = I±4£!>. (11) 

II qgA* Z 

Let 

D(x,y) = X(x)5 Xj y - J(x,y), (12) 

where 

y 

and suppose that we limit ourself to the ferromagnetic situation (e.g. 
J^{ x iV) > \J{ x iV)\i see [11)) expressed by the condition that the matrix 
D (|T2| ) is positive definite: 

D(x,y)>0. 

This implies 

D(q) = £ D(z, 0)e"^ = J 3 (0) - J(q) < £ D(z, 0) = D(0). 

z z 

Hence 

h > £[J 3 M) - J(z,0)] = D(0) > D(q) 

z 

and 

h - u(a 3 )[J 3 (0) - J(q)} = h- u(a 3 )D(q) >h- D(q) >h- D(0) > 0. 
From this, and from (|H]), it follows that 



0<u(F + (q)F-(q) 



Ua 3 ) < 1 



e 2/3{h-o;(o- 3 )D( g )} _ I — e 2/3h _ ]_ ' 
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The first inequality yields — 1 < to (a 3 ) < 0. 
Using (|TT|), one gets 

^ 2 - 1 + e^T[ < 13 > 

* + = + <") 

establishing a bound on the magnetization for low temperatures as a function 
of the interaction constants. The bound measures the deviation of the mag- 
netization from its ground state value, equal to —1, for small temperatures. 

Remark that for u(a 3 ) = —1, the ground state (/? — > oo) value of the fer- 
romagnetic system, the magnon creation and annihilation operators form a 
bosonic pair satisfying the canonical commutation relations 

and that the ground state a; is a magnon Fock state 

u (T + { q )T-{ q )) = hm emD{q l h] _ i = 0. 



In general, for all temperatures, the magnetization uj(<j 3 ) plays the role of 
the quantization parameter (a Planck's constant) (see (|7])) for the field of 
magnons. All quantum character of the magnons vanishes if one chooses 
the magnetic field h and/or the interaction constants (D(q)) and/or the 
temperature such that the magnetization vanishes. 

Concerning the magnetization fluctuation operators Fg(q), the ferromagnetic 
conditions (h — D(0) > 0, D > 0) are such that its infinite S'-limit F 3 (q) 
commutes with all other magnon observables (see (|])). They become clas- 
sical observables, and disappear completely from the action of the system 
Hamiltonian. It does not mean that there are no magnetization fluctuations. 

The original system with Hamiltonian iff in terms of the fluctuation ob- 
servables {ffiq), Fg(q)} qe h*, becomes in the infinite S'-limit a system of 
non-interacting magnons with Hamiltonian 

H A = £ t{q)F + {q)F-{q)i (15) 

qeA* 

where the magnon fluctuation creation and annihilation operators satisfy 

[T-(q)^ + (q')] = -uj(a 3 )5 q>q/ , 
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and where the spectrum is given by 



e(q) = 2 J 3 (0) - J(q) + 



Remark that if min(h, D(0)) > 0, then already e(q = 0) > 0, i.e. there is no 
condensation of magnons in the zero mode q — 0. In particular, this is the 
case under our assumptions. 

By inverse Fourier transform, i.e. 

II qeA* 

fll5|) can be written in configuration space: 

2h 

H A = 2^ D(x, y)F + (x)F-(y) + £ ^ (x)^(x), (16) 

with 

[^-(x),^ + ( 2 /)] = -a;(a 3 )4,, (17) 

and D is as defined above. Written this way, it is clear that 7^a is nothing 
but the 'quadratic boson Hamiltonian' of ||. 

Finally we look for the dynamics of the magnon excitation number operator 
J 7+ (x)J 7 ~(x). On the basis of (|TI|) the expectation value of this operator 
is related to the magnetization. Therefore the time evolution of this num- 
ber operator is related to the time evolution of the magnetic moment in a 
ferromagnet. 

The equation of motion for J r+ {x)J r ~{x) is given by: 
^T+{x)F-{x) 

= «[H Al f + (x)r(x)] 

= -2^(a 3 )£ J(x,y) (F + (y)f-(x)-F-(y)F + (x)) . (18) 

y 

Indeed, using (|P7D and J(x, x) = yields straightforwardly the result. 

This dynamical equation fll8|) can be compared with the macroscopic equa- 
tion of motion for the magnetic moment, first derived in |T^] (see also fl"3| , |I4| 1). 
A recent rigorous derivation of this equation for zero temperature is given 
T5J, applying a hydrodynamical limit or a Lebowitz - Penrose approxima- 
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tion ||16|| . Remark that equation fll8"|) however, is valid for all temperature 
equilibrium states. 
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